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ABSTRACT 


In this thesis we consider the problem of a charged, bound 
particle interacting with an electromagnetic field. 

First, we consider a classical model for the interaction of 
a bound, charged particle with an electromagnetic field. We determine 
the frequency shift of the oscillating particle due to the 
electromagnetic. field. 

Second, we consider a quantum mechanical model for a spinless 
bound particle interacting with an electromagnetic field. We determine 
the level shift, line breadth, and the probability of transition for this 
model. 

Finally, we apply the results obtained from the quantum 
mechanical model to the electromagnetic transition of an atomic and a 
nuclear system. From this application we form some general conclusions 
as to the dependence of line breadth and level shift on the binding 


potential. 
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CHAPTER I. INTRODUCTION 


If a system which consists of a bound, charged particle is 
placed in an electromagnetic field there will be an interaction between 
the particle and the field. Due to this interaction, the particle can do 
one of two things; first, the particle may absorb a photon and be excited 
to a higher state or, second, if the particle is in an excited state it 
may be induced to emit a photon and make a transition to a lower state. 
Due to the interaction with the electromagnetic field, the energy level 
of the particle is shifted to a new level, also when a transition occurs 
the line width is not a sharp line but has a finite width due to the 


radiative damping effect. 


The problem of radiative transitions has been studied by many 
different people since 1916 when Einstein“ presented his derivation of 
his now famous coefficients of absorption, spontaneous emission and 


3 really got 


induced emission. It was not until 1926, however, that Dirac 
on the right track in tackling this problem. In his work, Dirac presented 
the first satisfactory theory of radiation, in which he derived the 
Einstein formulas. 

The work on the quantum theory of radiation has been advanced 


by many people and this work has resulted in the treatises of Fermi“ and 


Heitler,> wherein reference to the earlier work can be found. 


The common approach to the problem of radiative transitions is 
through the use of perturbation theory. © In perturbation theory the 


interaction term is treated as a small perturbation on the system 
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of the electromagnetic field plus the bound particle. While the 
perturbation approach yields relatively good results, it would be better 
if the problem could be solved exactly. Unfortunately, this is an 
unrealized dream for Physicists at present. There are several reasons 
why an exact solution is preferred over a perturbative one. First, there 
is the aesthetic pleasure of being able to solve the problem exactly. 
Second, and most important, is that to use perturbation theory implies 
an assumption of analyticity in the coupling constant. This may or may 
not be the case; it is as yet unresolved. Third, there are cases where 
the coupling constants are quite large, i.e. nuclear forces, and a 
correct perturbation treatment is impossible. Since the complete 
Hamiltonian for a charged particle in an electromagnetic field cannot 

be solved exactly, we shall not use the complete Hamiltonian in this 


work, but a model which can be solved exactly. 


In the second chapter we consider a classical.model of a 
bound, charged particle interacting with an electromagnetic field. The 
particle is assumed to be bound in a harmonic oscillator potential and 
oscillates with a frequency, v. The shift in this frequency is 
calculated and is shown to be of the order on where e is the charge of 
an electron. 

In the third chapter we consider a non-relativistic quantum 
mechanical model of a spinless bound particle interacting with an 
electromagnetic field. The particle is assumed to be bound in a central 
potential, which may be the coulomb-field of the nucleus in the case of 
an electron; or nuclear field in the case of a nucleon. We calculate 


the line breadth, level shift and the transition probability. The line 
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breadth and level shift are shown to be of order et, where e is the 
charge of an electron. The results we obtain agree with the results of 
Chan and Razavy/ in their work on the theory of natural line breadth. 

In the fourth chapter, we apply the results of the third 
chapter to an atomic and a nuclear system. We give a general discussion 
on the dependence of the line breadth and level shift on the form of the 
binding potential. We find that in general a tightly-bound particle has 
its energy level shifted by a very small amount and its natural line 
width is quite small, while for a loosely-bound particle we find that 


the level shift and natural line width is quite large. 
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CHAPTER II. THE CLASSICAL MODEL 


In this chapter we shall discuss the classical problem of a 
charged particle, bound by a harmonic oscillator potential, interacting 


with an electromagnetic field. The Hamiltonian we shall use is® 
2 
H = (1/2m)(B(t)-efA(%,t)p(H)d2x]~ + (1/2)mv2R2(t) 
+ (1/8n) fitvxACZ,t)]* + 1612822, £) }a3x (2-1) 


caer ' ; P ; ; 
where p(r) is the charge distribution function which we take as 


o(f) = (62/4rr)e 8* (Z=2) 
where 
peel / eG (2-3) 


and where to is the mean square radius of the charge distribution. ro 
is assumed to be very small, approaching but never reaching zero, i.e. 
B is very large. ~In 2=i, A(z,t) is the vector potential, $(r,t) is the 
momentum conjugate to A(z,t), m is the mass of the particle, v is the 
frequency of oscillation of the particle, R(t) is the position of the 
particle, and P(t) is the momentum conjugate to R(t). 

The first term in 2-1 is the mechanical momentum of the system 
which contains both the energy of the particle plus the energy of 
interaction. The second term is the potential by which the particle is 
bound. The third term is the energy of the electromagnetic field. 


The equations of motion obtained from 2-1 are 
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a,A(z,t) = 6H/dd = 4nd(¥,t) (2-4) 
a,$(r,t) = -6H/6A = -(1/4n) curl curl A(Z,t) 

+(e/m)o(r) [B(t)-e fA(r,t)p (2) dx] (2-5) 
a,B(t) = -6H/éR = -mv2R(t) (2-6) 
d¢R(t) = SH/6B = (1/m)[B(t)-efA(Z,t)o(#)a3r]. (2-7) 


From 2-4 we obtain 
a2 A@,t) = 40 0, 8(,t). (2-8) 
Substituting 2-5 into 2-8 we obtain 


a2, A(z,t) = -curl curl A(r,t) + (4ne/m)p (t) [B(t)-efA(Z,t)o (Z)d?r]. 


(2-9) 
Using the coulomb gauge, where div A(z,t) = 0, we then have 
a2 AG,t) = V2ACE,t) + (4ne/m) (x) [P(t)-efA(t,t) 0 (F)d3x]. (2-10) 
From 2-7 we obtain 
92 R(t) = (1/m)3,P(t)-(e/m) fa, AG, t)o(r) dr. (2-11) 
Substituting 2-6 into 2-11 we then have 
a2 R(t) = -v?R(t)-(e/m) fa, A(z, t)o (r) 43x. (12) 


Equations 2-10 and 2-12 are the equations we now wish to 


consider. Let 


R(t) = [Re aie 
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>> : ‘ 
A(r,t) = JA er, (2-14) 
mn) 

Substituting 2-13 and 2-14 into 2-10 and 2-12 we obtain 
> 

(w*-vA)R, = (iwe/m) {A (Ep (F) d3x (2-15) 
> . 

(Vtw2)A CE) = -4neiup (ZR. (2-16) 


Since o(r), 2-2, is zero except very close to the origin we 
need to substitute in place of A) the average value of A (®) as r>0. 


We calculate the average value of A @) as follows. Let A) be 


given by 

74 r ' 

A (fr) = €,(w)B {sin w(L-r) }/r (2217) 
r 

where the sum is over the directions of polarization, and w = lk}, 


where k is the momentum of the field. We can write A () as 


Ke) = Bia) + Ai) (2-18) 
Feo : P Le 

where A (2) is the longitudinal part and A(z) the transverse part of 

Coe where curl Rd) = div A(t) = 0. The singularity of A(t) near 


meet A = = A 
the origin is 1/r therefore near the origin we can write A) = 


e, (w) /r; then we require div {8, (w) /r} = div A'(r). Since curl Re) = 0 
we can write KY) = vf(r); therefore 

div {é, (w)/r} = V7£(4) 

or 

£(¥) = (1/2)e,(w)°Vr. (2-19) 


From 2-19 we have 
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‘14 (1/2) V{2, (w)+V}r. (2-20) 


From 2-18 and 2-20 we obtain 


AN) = A@)-E@) = (8, (w) /}-(1/2) 06, (w) +e 


it 


(€,(w)/2r} + {[€, (wu) +] /ar3}z. 
Consider 
2, (w) ACE) = (1/2r) (1 + (8 (w)-F}2/r2]= (1/2r) [1 + cos26] (2-21) 


where cos@ = & (w)*t/r. Using 2-21, the average value of @) (w) AT (x) is 


then 

ane a 

e, (uw) -A' (x) = (1/2r) J" (1 + cos?6]sindd6/{"sinede = 2/3r. (2-22) 
0 


We therefore write 


> > mn > 
A (x) = )(2B /3r)e, (w)sinu(L-r). (2-23) 
W i W 


Substituting 2-23 into 2-15 we obtain 


(u2-v2)R, = (2ive/3m) JB’, (w) fo (E) {sinw(L-r) /r}a3x (2-24) 
A 
or 
R = (2iwe/3m) (1/w2-v2) )B*é, (w) fp (£){sinw(L-r) /r}d3r. (2-25) 
U) , wr 


Substituting 2-25 into 2-16 we obtain 


> > 
(V2tu2)A (x) = {81w2e2p (£) /3m2 (w2-v2) })B*€, (w) fo () {sinw(L-r) /r}d3r. 
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(2-26) 
Substituting 2-17 into 2-26 we obtain 
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Substituting 

V*{sinw(L-r)/r} = -w2{sinw(L-r)/r}+sinw(L-r)V2(1/r) (2-28) 

into 2-27 we obtain 

sinw(L-r)V*(1/r) = {8ne*w?/[3m(w2-v2) J }p(£) fo (%) {sinw(L-r) /r}d3r. (2-29) 
The integral occuring in 2-29 can be written as 


fo(#) {sinw(L-r)/r}d3r 
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sinwL{B/[1+(w2/82)]} - weoswL{1/[1+(w2/87) ] 


(2-30) 
where we have used 2-2. 


With the assumption that 8 is very large we can neglect the 
second term in the denominator of each of the terms in 2-30, we then 


obtain 

fo (©) {sinw(L-r)/r}d3r = {sinwL/r)} - weoswL (2-31) 
where we have used 2-3. Since V*(1/r) is given ee 

v2(1/r) = -4n8 (x) (2-32) 
we can write 2-29 as 

-sinw(L-r)6(t) = {2e2w2/[3m(w2-v2) ]} { (sinwL/r,)-weoswh]o (x) (2253) 


where we have used 2-31. Integrating 2-33 over all space we obtain 
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-sinwL = { 2e*w* /[3m(w2-v2) ]} { (sinwL/rg) - weoswL] (2-34) 

where we have used the fact 

(82/4n)f(e 8" /r)d3x = 1. (2-35) 
From 2-34 we obtain the eigenvalue equation for w, i.e. 

tanwL-2e*w? /{ 3m(w*-v?)+2e2w*/ry} = F(w) = 0. (2-36) 


We note that if e*0 in 2-36 we obtain the eigenvalue equation for the 


uncoupled system, i.e. 
Veeeant/ism = 05152)... (2-37) 


We now want to show that for each V there is one and only one 
W. Consider 2-36, where nn<wL<(ntl)7. There are two cases to consider 


in studying this equation. First there is the case where n is such that 
3m(wé-v*) + (e7uc/r,) = 0, (2-38) 


i.e. the resonance condition. Second we have to study the case for all 
other values of n. First we consider the second case. If we draw a 


sketch of F(w) versus w, where 
F(w) = tanwL-2e7w?/{3m(w*-v?)+2e*w*/r} 


For nr < oL = (ntl) we obtain 
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Figure 1. 


where 


. 2e?[nr/L]° 
3m[ (nt/L)2-v2] + (2e7/r))[nn/L]? 


= 2e2[ (nt1) 1/L)3 


Beale i/o SINE Gi oe, UE aN hy yin 

3m{[(n+1)a/L]*-v~} + Qe*/r)[ (ntl) 1/L] 
We see from figure 1 that 2-36 only has one root in each interval 
nt < wl < (ntl)7. We also note that the root of 2-36 satisfies 


nt < wih < nn + 7/2. We also notice that as n > ~ the value of 


ot/L + n/2L. 
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We see from 2-38 that for resonance Wp satisfies 
ily 
bp =O) [ler 4 2e"/ Snr) ] 


For this case we obtain the following figure 


Figure 2. 


iZ 
From 2-36 for each Wn such that 2-38 is not satisfied we 


obtain 
tanw L = 2e*w3 /[3m(w2-v) + 2ertie/ iad: (2-39) 


It is important to note that v in 2-39 is not one of the VV, a0 2437. ¥ 
is the frequency of oscillation of the particle due to the harmonic 
force in the absence of the electromagnetic field. The VyS are the 
frequencies of the electromagnetic field in the absence of coupling to 


the particle system. The wis are the new frequencies of the coupled 


n 
system, w) is the shifted frequency of the oscillating particle and the 


other w 8 are the shifted frequencies of the electromagnetic field. 


From 2-37 we have 
tanvL = Q. (2-40) 


From 2-39 and 2-40 we then obtain 


tan(w-v,)L = 2e7w3[3m(w2-v*) + 2e*w2/ro] (2-41) 
or 
wmv, = (1/L)Arctan{2e*up/[3m(wp-v*) + 2e*wn/ro]}- (2-42) 


Since the argument of Arctan, in 2-47, is very small we can write 
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(1/L) {2e we /[3m(o, v“) + 2e we/roli. (2-43) 
For n = O we obtain from 2-43 


Cage (1/L) {2e%w3 /[3m(u§-v*) + 2e7wi/r,]}. (2-44) 
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Since we are mainly interested in the shift in the oscillating particle 
frequency we will not consider the equation for n = LEZHS Bs he 


From 2-44 we can write 
wg =v? - [2e¢we/(3mrg)] + [2e2w2/(3mL) ] (2-45) 


or putting Wy = Vs since Wy-v % e*, on the right hand side we have 


Wo = vite[2e-7 (nr) ] + ree (amnyn (2-46) 


Since ry is assumed to be very small we can neglect the last term in the 


parentheses in 2-46, we then have 
Wo = v[1-2e2/(3mr,)}2/? (2-47) 


or if Ee 2e2/3m we can write 


ip = =le-/ (ne) 10 (2-48) 


We see that the frequency shift is of order e* and proportional to v. 
We notice that if r, is decreased to zero in 2-47 we obtain 

imaginary values for Wo « These imaginary values correspond to the 

phenomenon of Jel eeceaieratiody.. We see that the minimum value that 


Lo can assume such that Wo is real is 


ry = 2e2/3m. (2-49) 


The classical radius of the electron is e*/m, thus we see that to is 
approximately equal to the classical radius of the electron. 
The phenomenon of self-acceleration arises because there are 
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first considered by Dirac, which eliminates these unphysical solutions 
is that the acceleration goes to zero as time approaches infinity. 

It should also be noted that no classical theory can explain 
the structure of the electron and therefore all structural concepts, 
such as r@dius, should not occur in the classical theory. The 
structural properties should be explained, if indeed they can be 


explained at all, within the framework of quantum ene 
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CHAPTER III. THE QUANTUM MECHANICAL MODEL 


In this chapter we consider a quantum mechanical model for the 
non-relativistic interaction between a bound, spinless, charged particle 
and the electromagnetic field. An example of such a system is the 
nucleons in the nucleus interacting with an electromagnetic field. We 
shall consider only processes involving a single charged particle with 
the emission or absorption of one photon. 

In section 1 we discuss the model and obtain the Hamiltonian 
we shall be using. In section 2 we derive the equations of motion 
which we solve in section 3 and obtain expressions for the level shift, 
line breadth, and the cross-section for emission or absorption of a 


photon. 


1. The Hamiltonian 
The complete Hamiltonian for a non-relativistic, bound, 
spinless, charged particle interacting with an electromagnetic field 


is,* (h=c=1) 


Hoe 32 Hi S-) 
where 
Hy = [{(1/2m)VyeCE,t) Ww, t) + V(E)Va,t)y@,t) }dex 

+(1/8n) [{B2(z,t) + [vxACr,t) ]?}d3x (3-2) 
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Hy = (de/2m) [Ar ,t) + (yee, t) H(z, t)-Wy*(z, t) H(z, t) Jase 
+ (e2/2m) f[A2(z,t) va(z,t)v(e,t)d3r 


+ (e2/2) ffyr(e,t)wr(ese)y (est) v(z,t) /|r-2| d3ra3r-, (3-3) 


The first term in 3-2 describes the charged particle, mass m, 
bound in a potential V(r). The second term in 3-2, where we use 
Gaussian units, is the radiation energy of the electromagnetic field. 
Both E(z,t) and A(z,t) are transverse fields. 

The first and second term in 3-3 describe the interaction 
between the charged particle and the electromagnetic field. The third 
term in 3-3 is the electrostatic self-energy and can be absorbed into 
V(r) in 3-2. 

In the model we shall study we modify the first and second 
terms of 3-3 so as to obtain a Hamiltonian which can be solved exactly. 
The justification for these modifications is mainly that the resulting 
problem can be solved exactly. Another justification is that this 
model explains some of the experimentally observed properties of the 
interaction of bound charged particles and an electromagnetic field. 
Also, once we have obtained an exact solution of our model it would be 
possible, not done here, to apply perturbation theory and thereby solve 
the complete problem, where the perturbation used is the difference 
between our model and the exact Hamiltonian. 

The Hamiltonian in 3-1 is the classical Hamiltonian, therefore 
it must be symmetrized before it can be used in the Quantum Mechanical 
formulation of the problem. This symmetrization is needed since the 
Hamiltonian, or any observable, is required to be Hermitian. We will 


therefore have to symmetrize the Hamiltonian we use for our model. 
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The Hamiltonian we shall consider is 


7s ae (3-4) 
where 
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(3-6) 

where we have symmetrized H so that it is Hermitian. In 3-5 and 3-6, 
i and f refer to the initial and final states of the system. AY ,t) 
and A (f,t) are the positive and negative frequency parts of the 


electromagnetic field. Assuming periodic boundary conditions and unit 


normalization volume we can write,l 
| le) 
> k,A > > 
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kd ee Pree 
where the photon energy is k = | & | and é, (k) is the unit polarization 


vector orthogonal to k. The sum over k extends over all plane wave 
momentum states k in the unit normalized volume. The sum over A extends 
over the two allowed directions of polarization of a photon having a 
momentum k. ay (t) La, ,(t)] is the creation [annihilation] operator 


ul 
for a photon of momentum kK and polarization 8, (k). We also have, 
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> > 
A(r,t) = Att) + FG, t) 


(3-8) 
E(z,t) = - BAC, t)s 
Defining ULC) as the eigenfunction of the Schrédinger 
equation 
((-1/2m) Vv? + V(z)}U,(£) = EU, (2) ee 3 (3-9) 
we then can expand v, (Ft) and y*(r,t) as 
a 
> am > 
y (f.t) = bi (t)U, Ce) 
(3-10) 
k(r are 16 
ye(z,t) = bo (t)UX(E) 


where bi (t) [b (t)] is the creation [annihilation] operator for a 
charged particle in the eigenstate a. 
The b(t) and ay , e) satisfy the following commutation 
3 


relations 
[b (t),b'-(t)] = 6 
Qa Qa + a 
(3-11) 
la, Ct) sayy -(t)] = Jal ae ee 


with all other pairs commuting. 


Substituting 3-7 and 3-10 into 3-4 we obtain the following 


operator form of our Hamiltonian 


or 
. ia 


(@£) tt =< » Cig = Cogulan 9 


en (a. 3)49 bre (3.7) basqxs aso modd 8 


| i 

) UG) d = Ge 

a 

(OL-£) i 

+, ’. Gea 

(A)auG) a = — * 

‘vin 
o 


6 102 totaxesqo [acizslidtaos] notisexo sda el (3), dj ieee rar 


4g) 
» s3etensgts sda nt sfoisisq beg 


netsasummoo goiwolfot ata yierise (3), 48 bas (3) 3d oaT 


‘ays | dob \, 
~ca? 7 SO ee 


oe 
(11-8) 
<7 soa, 90 = TOD — r 
: 
, 


eter ns 


19 


eet 
H = Epb (t)by(t) + B,bi(e)b4(t) + Di Lay (e)ay ()+(2/2) 
+ e) (1/2k)1/? igceypt(eybe(ta (t)+8*(k)al (t)b/(t)b 
L . F « t)b,-(t)b, (t)] 


+ (me2/2m) YJ (1/kqy* ToT Ceye (eyB, Cea) 2ah (eda, (£)48 4] 


q,k a=i,f 
(3-12) 
where the sum over } is suppressed, and where 
(je =e Gin 4 de a > >\> >, 43 
8 (i/m) (47) fexp(ik-£)U, (F)é(k) -W,(4)d r (3-13) 
and 
B (ka) = &(k) +8 (q) fexpli(K-q) -r]U*()U, (r) der. (3-14) 


If the last term in 3-12 is dropped we obtain the same 


17 in connection with the 


Hamiltonian as studied by Chan and Razavy 
problem of natural line width. 
The Hamiltonian 3-12 is in the Heisenberg representation; we 


can transform to the Schrédinger representation by the transformation 


of operators given by 


O(t) = exp(iHt)0(t=0)exp(-iHt). 


Thus in the Schrédinger representation the Hamiltonian has the same 
form but with all operators evaluated at t = 0. We shall drop the t 
parameter and assume all operators are evaluated at t =0 and use the 
time-dependent Schrodinger equation to study the equations of motion in 


the next section. 
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2. The Equations of Motion 
We shall assume that the initial condition of the system is 


given by 
|y(O)> = |0,1,> (3-15) 


i.e., the charged particle is in its initial state, with no photons 
present. Due to the form of the Hamiltonian, 3-12, |y(t)>, the state 


after time t, must have the agen” 


ie} eo= oC, (6) )0.15> 4:4) 6. Ge) [dp dee (3-16) 
k 


where [Lae denotes the state with the charged particle in its final 
state and one photon of momentum K and polarization é, (k) present. 


Comparing 3-15 and 3-16 we see that at t = 0 we have 


Cy Ct =0) =1 C(t = 0) = 0. (3-17) 


1S) 


Following the method of Haake and Weidlich”’~ we write the 


Schrédinger equation 
ia ,|y(t)> = H] p(t)> (o=13)) 


where H is given by 3-12 but with all operators evaluated at t = 0. 


Substituting 3-16 into 3-18 we obtain 


da WY LpZ 
fagGp(t) = ByC)(t) + JLeB ()/(20 ellos (3-19) 
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EERO | 
9, (t) = (Ep + k)C,.(t) + [ep*(k)/(2k) JC, (t) 


+ (ne2/m)) (1/ (ka) */7 1B (-W) 6, (€) (3-20) 
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where 

Ey = E, + (ne2/2m) BS) (3-21) 
and 

ee (me*/2m) )[B -(0) /k] + (1/2))k. (3-22) 


k k 


The two divergent quantities, 3-21 and 3-22, are the 
renormalized energies of the initial and final states of the system. 
The expression 3-19 and 3-20 are the equations of motion of our system. 

To find the normal modes of this coupled system of 


j ; 2 
equations we write the Fourier series for Cy (t) and C(t) as 2 


noe Lp Ca) g, (ad expl i (a ESE) n = 0,k (3-93) 
where 
pia) = C¥(a)C,(t = 0) + pepe: = 0). (3-24) 


Substituting 3-17 into 3-24 we have 
p(w) = CX(w). (3°25) 
Substituting 3-23 into 3-19 and 3-20 we obtain 


yack) C, (w) (3-26) 


we, (w) errytan) 
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1/2, 


wO,(w) = (By + k - E,)C (w) + LeB*(K) /(2k)°! IC, (w) 
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By eliminating C, (w) between 3-26 and 3-27 we obtain an 


equation for C, (w) 


C, (w) = ee ett) (3-23) 


where 
K(g,k) = 1/[ CE g+k~Ey—w) (a) */?{ (8% (k) 8 (2) 1 /2u#(n /m)B.(2-K) } (3-29) 


In the next section we consider 3-28 and its solution and 


determine the level shift, line breadth, and cross-section. 


3. Level Shift, Line Breadth and Probability of Decay 


We can write 3-28 as 
C(k,w) = ~e? fO(L,u)K(L,k5w) at (3-30) 
where we have used the relation 
d + [2/(2n)3]faek (3-31) 


which is valid when the normalization volume tends to infinity, the 


factor of 2 in the numerator is due to the sum over the polarization and 


where 


re u 2 1/2, ,B*(k)B(2) 4 7 e 
K(2,k3w) = 22 (a LS ey) [Re (Cee aad om + -B,(2-k)]. 


(3-32) 


We now write 3-30 as 
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where we have used 3-32 and where 


=e (fu) foca, w)B CL) 2 os dy. (3-34) 


Substituting 3-33 into 3-34 we obtain 


A= [ f°f" Clasw)Be(a-e) a (ua? dads 7, e2 f*|a(ey|2eazy72. (3-35) 


Og hee Ey - 4) out 2n2 OE gtL-Ey-w | 


We now define the reciprocal kernal, T(2,k3;w), for 3-33 such that 


T(g,k3w) = K(g,k3w) - e°J T(L,45u)K(q,k3w)dq (3-36) 
where 

Be (2-k) 24 
K(2,k;w) = ci ee cee (3-37) 


rm (2k) 1/2 Cp +k ~ E, - w) 


We then obtain from 3-33 and the reciprocal kernal, T(2,k3w), 


the following 


Zax po 2ax(L A 
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(3-38) 
From 3-34 and 3-38 we obtain the eigenvalue equation 
y= - 22 pla lfade yet pp Bue aan w) 9?! Faqs (3-39) 
272 9 E,th-Ey-w ogc (q) 1/2 (H Egtq-E,-w) 


Although, in principle, we can solve 3-38 and 3-39 and obtain 
the desired results; i.e., level shift, etc.; we shall first consider a 


simplified version of the problem which can be solved very easily. We 
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24 
can then perhaps apply the knowledge gained in solving the simpler 
problem to the one above. 

Our simplification consists in setting BC e-k) equal to zero. 
This is only justified by the fact that the resulting equations can be 
solved very easily and exactly. We then have, instead of 3-26 and 3-27 


the following 


WG (w) e) 11/28) "718006, (0) 


(3-40) 


we, (w) = (Ey +k- E7)C, (w) + [ep*(k)/(2k)/ 7}, (w) 


where now instead of 3-21 and 3-22 we have 


foi las Seu 


We also have, using 3-30 and 3-32, 


co 


C(k ,w) = -e*f AGES EI SIGCIEEE ° (3-41) 


OF a2 . ee 
27 (E,tk EL w) (kq) W 
The eigenvalue equation 3-39 now becomes 


wo + (22/20) [| 8) |2/ Ch tk-B,-w) Idk = 0 (3-42) 


We now normalize Cy (w) and C, (w) by requiring 


Ch (w)Cy(w*) + LOL(w)C (w*) = 6, ye (3-43) 
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From 3-40 we find 


Z 
j ]}C, (w)- 


Zz 


£ 
Cw) = feg*(k)/[ (wa - BE, ~ k + EL) (2k) (3-44) 
Substituting 3-44 into 3-43 with w = w°, we obtain 
-1/2 
Cy (w) = {dA(w) /dw} (3-45) 
where 
2 ees 2 : 
Aw) = w+ 22 f_ |B) |Pkdke 3-46 
ee Creidn Seay eng a eae’ 
Substituting 3-45 into 3-44 we find 
aa fe) 
O.(w) = feB*(k)/[(w - E, - k + £,)(2k)4/2] [aa(w)/aw]/?. (3227) 
From 3-23, 3-25 and 3-45 we then obtain 
-1 
Co(t) = } exp[-i(w+E,)t] [dA(w)/dw] (3-48) 
w 
and 
yee awl Gass 
c(t) = ){eB*(k) /[ (2k) “ (w-E,-kt+E;)]} [dA(w)/dw] “exp[-i(wtE,)t]. 
Ww 
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Letting z = wtE -E we can then write 3-48 as 


Co(t) = [exp{-iE,t}/2ni] $ {exp(-izt)/A(z) }dz 
(6; 


(3-49) 


(3-50) 


where the contour C contains all the roots of A(z) and A(z) is defined 


by 3-46. We must now consider A(z) if we wish to evaluate 3-50. 
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Letting 


Z=wt ET = Ep (3-51) 


we can write 3-42 as 


Mzy =z = E+ Bt (eP/2n?) JE] 8¢k) |?/ (k - 2) ]kdk (3-52) 
Expressions 3-50 and 3-52 are the same equations as studied by Grants 
and we shall employ the same methods as Chan to solve them. 

We can analytically continue 3-52 onto the complex plane if 
|2 24 


we can analytically continue |8(k) We shall assume that |8(k) |? 


can be analytically continued and we then obtain 


a\(z)=2-E, +E, + (27/21?) f °C | 8k) |2/( ee (3-53) 


where 3-53 is defined on the first Riemann sheet for 0 < arg(z) < 27. 
It can be seen that 3-53 has a cut along the positive real Perea 
The integral in 3-50 can then be obtained by integrating along the 
contour, Figure 3, in the first Riemann sheets” 


We must find the zeros of At (z) in the lower half plane, 


excluding the cut. Substitution of z = x + iy into 3-53 yields 


(3-54) 


x = Ey ~ Ey - (22/20?) J [|B (k) |7/ Ck - x) ]kdk. 
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Figure 3. 


Since we expect 3-50 to behave as a damped exponential which 


arises from a complex root of A(z) with a negative imaginary part we see 


from 3-54 that we must look onto the second Riemann sheet for the 


desired root. 


Analytic continuation onto the second Riemann sheet, for 


0 > arg(z) > -2n is effected by defining for real z = x, 
at ex - ie) = At (x + ie). 
Then for z = x + iy we have 


att(z) = * (2) + (ie2/1)z|8(z) |. 


We can now evaluate 3-50 by using the following contour. 


(3-35) 


(3-56) 


€ siugtt 


dotdw Lalsasmogxs boqmsb' s 2s oveded o2 O@-€ toeqxe ow sonte | Pr 7 
398 aw, 38g qaantgemt ovitegom s ditw (s)d to Joor xslqmos s mort cet 


7) 
edd tot desde ansmatd broose sd3 ojo sink deum ow tarts sent at 
ao) | a 
is ado Bab a! ; 
toi" a anemsti baoase sd3 ojo notiswatinos ottylenA . Le. 


7 , - 4 


«x # s Iso1 103 gaiatieb yd bstastis ak 7S- < (s)gts < o 
: ; ee 7 


(22-€)_ (at + x) te = Gr = yy 
, ee 
7 TS ee. 

coe 3 ) oved ow yt + x = 8 10% aod 


a 
(22-£) | | eee + (sts = @ " 


4 Pin mi ss coca a ee 


28 


Figure 4. 


The solid portion of the contour is on the first Riemann sheet while 
the dotted portion is on the second. We now must find the zeros of 
Atl (z) in the second Riemann sheet. We assume, for simplicity that 


IT 2) has only one solution at 


Za = E, = ED Smecirlyy) 2)% (9-579 


Substitution of 3-57 into 3-56 and then equating the real and imaginary 


parts of 3-56 we find 


Boe + Sf 160) ke + B= By)dk ~ Son = 0 (3-58) 
TT 
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where 


Se te (i/2)]| BLE, = EBs ar 2yF (3-60) 


Equation 3-58 gives the level shift, E, - EL: and 3-59 gives the line 


R 


breadth, [3 these expressions are the same as the ones obtained by Chan 


and Ragavy .-- 


Approximate solutions for E 


ys EL and [, which will be useful 


later on, may be obtained by assuming ER ~ EL or [ > O we then obtain, 


E. - E. = - ef pf |B (k) | ®kdk (3-61) 


R I 210 E 
27 k + By, ED 


and 


2 
ee < (Ep - Ey) |B(E, haar (3-62) 


The results in 3-61 and 3-62 for E, - E,_ and I agree with that obtained 


R ale 
by first order time-dependent perturbation een era lepine The symbol P 
stands for the principle value of the integral. 
We can now write 3-50 as 
= exp(-iEpt) exp (-izt) = 
C,(t) aed : {f{ + f} Sarre Zz (3-63) 
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where C, and Cy are shown in figure 4 and A(z) is either A (z) or A (z) 


depending on whether the contour is on the first or second Riemann sheet. 


The integrals over Cy and Cy are then, 


exp(-iEpt) exp(-izt) exp[-iE,t-(I/2)t] (3-64) 
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and 


exp(-iE,t) ;exp(-izt) , = Oxp(-iEpt) , -oexp(- ixt) , 
2ni Cc, A(z) Qri wo Atl (x) 


ou 


rae “exp (- ixt) gy} 


il 
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(3-65) 
then we can write 3-50 as 
exp[-iEpt-(r/2)t 2 2 a 
c,(t) = = R ] z ee (Ae as xt) cux xp(-i8,t) 
id AS AZ) P21. a amaqectea: 272-0 AEGAN Gx) 
z=E,-E., Ciy:2) (3-b¢) 


where we have used 3-56 with z = x (real) for the integral over Coe 


Since we have set Be (g-k) equal to zero we can write from 3-20 


* 

a) ~LB2 CE) exp l-i(ktEy )t] if exp[i(ktE,)t]C,(t)dt (267) 
(2k) 

where we have used 3-17. We could now substitute 3-66 into 3-67 and 

solve for C(t), however, the solution is so complex that no useful 

conclusions can be drawn from it. Therefore to obtain a useful 


expression for our results we write C(t) as, 


exp[-iE t-(1/2)t] . 
Cy (t) ee : (3-68) 
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From 3-56 we can write 


i 20 
dA : 5 ie? 
ia) ee - yf 18) aS ees Cea (3-69) 
ae 272 0 (k-z)2 1 dz 


thus we see that a good approximation for ee)! ee 
Zz 


dall(z) 7 


tks = 
dz (3-70) 


Substituting 3-70 into 3-68 we then have 
Co(t) = exp[-iE,t-(I/2)t]. (3-71) 
Substituting 3-71 into 3-67 we obtain 


eg*(k) exp[-iEpt-( r/2)t]-exp[-i(k+E_)t] 
ee R E hs 
C(t) ae ( oi: (Sie) 


[E,, +k - ER + Cie 2) 
Now that we have found C(t), 3-72, we can obtain the cross- 
section and the oscillator strength. The oscillator strength is given 


ee 


a Dee a 2 = 
£ a(t) = (m/2ne*t) J |e, (e) | dk (3-73) 


where the integral in 3-73 is over the energy variable. From 3-72 we 


find 


2 2 - 1-2 St why 74 k+E_—E,,)t 
puieas: RCO esl ie lise EC a, 
k 2|k| Cee ey 


ive ee 
where we obtain [k| 5 in the denominator, from [(k) ]¢ since k can 


assume both positive and negative value. Defining z = k + Ep ER 
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and substituting into 3-73 and substituting the results into 3-74 we 


obtain 


ih = 2 
ae AUC ae 
4n2 t (2-EptEp) z*+(T/2)2 


{exp(-I't)+1-2exp[(-I'/2)t-izt]}dz 
(3-75) 


where the integral in 3-75 is done over the following contour. 


Figure 5. 


The contour of figure 5 is chosen in such a way that it includes only 
the singularity at z = (-il/2). The selection of this particular 
contour is justified since in the limit [ > 0 the results we obtain for 
aes see 3-81, is the same as that obtained from perturbation theory. 


We then obtain from 3-75 


=Ee=(@0/2)iy el 
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f Siena Ree FUEL) 
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a3 
We see that foe 3-76 is time-dependent; this time-dependence comes 
from the finiteness of I. Taking the limit [ > 0 we see that 3-76 
becomes time independent as is normally expected. 
We know from perturbation theory that frp = 1, we can 
obtain this result by assuming I > 0, Ep = EL» and by considering the 
electric dipole transistion; note that in this approximation E, = E, 


I al 
and E,, = Eg- We use 3-13 to calculate 8(k), 


B(k) > [iC4n)4/?/m} f up) a up@)a%x 
Zz 


Ub he 


= -[(4n)~/“/m] < I | P, | F > ; (Sad ) 


where De is the mechanical momentum of the charged particle. We have?” 


1 | P, es (Ey Eee acne (3-78) 


Substituting 3-77 and 3-78 into 3-75 and taking the limit as I > 0 we 


find the following 
Lo = ) 2m (E, eh), |< I | Zz |F ae (3-79) 
L I 
or since we have the following, 

; alors = 
une E,) |= 2 | 2 [rele = 1 (3-80) 
I 
we obtain the desired results, 
; = (3-81) 
cee i 
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for one charged particle. 


One finds that if fir denotes the time average of fet over 


the time interval 0 < t < (1/T) then 
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fip = (-1) j Se Pea0)) = 3/905 fowl 0) (3-82) 
ee ee Le IF 
where 
C1yT) 
fig =1 Y, Fi p(tdt (3-83) 
and 
alr) : na epee 
s {l-exp(Tt) }t ode =(-1) ) oo) (3-84) 
n=1 . 


and fi p60) represents the oscillator strength evaluated at t = 0. 

We can now calculate the cross-section for the emission of a 
photon of momentum k. The matrix element for the photon to be emitted 
with momentum between k and ktdk is given in the Schrédinger picture by, 
M, p(k) dk = kPdk < dp(k,t = 0) | H(t = 0) | oz Ct = 0) > (3-85) 

Tr 


where 


ese Ce) | ta 


i 
(3-56) 
op (e)> 965 ¢t) | 0,1, > 
and 
H“(t=0) = e) (2/24) 1/718 (bib ga 48*CH) 3.054) (3-87) 
Substituting 3-86 and 3-87 into 3-85 we obtain 
(3-88) 
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Substituting 3-71 and 3-72 for C,(0) and CE (eo) we find 


2 4 252 
[Mp p(k) |2 = a ee EUS OF Oe (3-89) 


Ar (Ep-Egtk)2+(r /2)? 


The flux of photons per unit frequency interval per steradian 


for a given polarization is 2/202 ,>> thus for both polarization 


directions the flux is k2/n2. The cross-section can now be determined 


by the relation, 
ek) 2m tk) [=Cie 7k). (3-90) 
Substituting 3-89 into 3-90 we obtain 


4 
meee ee CUCU ea (3-91) 


2m (Ep-E,tk)*+(T/2)? 
Assuming 
Bo = Bo = Eo = Be 2" (3-92) 
we find, using 3-62, 
T = (e2/7)k| BC) | 2 (3-93) 
then substituting 3-92 and 3-93 into 3-90 we obtain 


Tr T2 , 
pot a eee ees Care Lee ee (3-94) 
a) oe ya (Ey-B, tk) ?+(1/2)? 


: ; 36 
Eqauticn 3-94 is the well-known resonance scattering cross-section. 


The factor of 2 in the denominator arises from the two polarization 


states of the photon. 
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We now must go back and consider 3-39 and attempt to determine 
the level shift and line breadth as given by this equation. We shall 
again assume that there is a single solution of 3-39 on the second 


Riemann sheet and that this, as before, is located at 

z= ER7E,- GP /2) (3-95) 

or by using 3-51 

wo = EE, (it /2). (3-96) 
Analytically continuing 3-39 onto the second Riemann sheet, 


and substituting 3-96 into the results, we obtain 


ez ms B(2) |22de 


70 a fe 
20 OE te Ett /2) 
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| nee oo 3/2 
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(3-97) 


Equating the real and imaginary parts of 3-97 we obtain 


eee ee es eee 
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where 
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(3-100) 


These results for the level shift and line breadth, 
3-98 and 3-99, differ from the previous results, 3-58 and 3-59, obtained 
with a simpler model, by a term proportional to eae 

Expressions 3-98 and 3-99 along with 3-38 are the solution 


to our problem. While these results are not simple, they can be solved 


by numerical methods. 
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CHAPTER IV. DEPENDENCE OF THE LEVEL SHIFT AND 
LINE BREADTH ON THE BINDING POTENTIAL 


In this chapter we shall determine what, if any, effect the 
shape of the binding potential has on the level shift and line breadth. 
We shall use for the level shift, AE, and line breadth, I, the 
expression 3-61 and 3-62. 

We shall only consider the case where both the initial and 
final states are bound. 

The dependence of the level shift and line breadth on the 
potential is through the dependence of g(k), 3-13, on the wave functions 
U, (F) and U_(t), which in turn depend on the potential. In the first 
section we derive expressions for the level shift and line breadth for 
general forms of U,(¥) and U,(#). In the second and third sections we 
have assumed specific form for U, (fF) and Us (¥) and from these determined 


the level shift and line breadth. 


1. Calculation of 8(k) and General Expressions for AE and I. 


From 3-13 we have 
B(k) = (i/m) (4n)/? fexp(ik-£)U, (x)e(k) + WU, (r)d5r. (4-1) 


Taking y as our polar angle and taking e(k) in the z direction, i.e. 
eth) =a(0.0, 1.) and k in the x direction, i.e. k = (k,0,0) we then 


obtain 


eo = (i /m) (47) */? fexp (ikrsinécosy)U, (F) cosé (3 /dr)Ug (F) sindr*dydedr. 


(4-2) 
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For the sake of simplicity we will consider transitions from a P state 


tO an: S state, i.e. U, (f) and Us (£) will be written as 


U, () = (3/41) 2 cosen,R, (x) 
U,() = (1/4)+/?neR, (x) 
where 

ae “1 
Nf = ae [Pee enade 


NE oF [Re (x) |2x2ar] 7 


Substituting 4-3 into 4-2 we then obtain 


B(k) = (4 /m) (3/40) 5 Ne fexp (Akreindcosy)R, (x) [(8/8r) Re (r)] 


cos*Osinér2dydédr. 


The integral over y is given by?! 


27 
{ exp(ikrsin§cosy)dy = 27J,(krsin6). 
0 


The integral over 6 is then given by 


AL i 1/2 
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We now can write 4-5 as 
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This is the desired expression for B(k). 
From 3-61 and 3-62 we have 


co 


AR Sen, "= Re =oeet pp: 160c) | 2kdk 


207 "0  k+ER-E, ee 
Tr = (e*/n)(Ep - Ey) |B(Ep - E,) |? (4-10) 


where from 4-8 we have 
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(4-11) 
From 4-9, 4-10 and 4-11, with an assumed form for R, (x) and 
R(x), the expressions for AE and [ can now be calculated. We now turn 


to two specific examples of this calculation. 


2. The Atomic System. 


We shall assume that the normalized radial wave functions 


R, (r) and Re (4) are given by 


R(x) = {20°/2/(3) 4/7] rexp(-nx) 
(4-12) 


R, (4) = 2y3/2exp(-vr) 


where » and v are proportional to the energy of the particle and are 


given by?” 
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y= (2E,/e*) 


We notice that the dimensions of uw and v are fm if energy is measured 


in Mev and e* has dimensions Mev fm. 


Substituting 4-12 into 4-8 we obtain 


B(k) = eat GG yyl2r “312° Sexpl- (utv)r]Jdr. (4-13) 
m 0 (kr) af 2 


The integral over r in 4-13 is given by? 
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Substituting 4-14 into 4-13 we obtain the final form of 8(k) 
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Substituting 4-15 into 4-9 we obtain 
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Substitution of 4-15 into 4-10 gives 
2 
T = pos (uv)? Wistesn (4-17) 
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Putting H and c in 4-16 and 4-17 where they should occur we obtain 
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256Ke* e (Ep~Ep/Uc) 
= (wv)§ ———4 = —___, (4-19) 
c“m { (uty) *+[ (E,-Ep)*/K2c7]} 
then AE is given in units of Mev and TI in units of sect, 
We shall assume that ER = Ey since the difference ER-Ey Se Or 


order e? and in this approximation correction of this order can be 
neglected. We notice that 4-18 and 4-19 are symmetric in y and y; 
therefore we shall only consider the case where the final energy, i.e. 
v, is kept constant and the initial energy, i.e. uy, is allowed to vary. 
The fact that 4-18 and 4-19 are symmetric is not in general true, as 
will be seen in section 3. It is obvious, from 4-8, that 8(k) is not 
in general symmetric; thus 4-18 and 4-19 are not in general symmetric. 


Considering 4-18 and 4-19 as u > ~ we find 


AE v (1/y') 

(4-20) 
To (1/5) 
For u > 0 we find 
AE v yp? 

(4-21) 
rv The 


Since both AE and I are finite everywhere we then obtain the following 


figures. 
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Figure 7. 
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Therefore, we see that if up is large, i.e. a strongly bound 
state, thenAE and [ both are small. This is due to the fact that the 
particle is bound very tightly, hence the effect of the electromagnetic 
field on its energy is very slight. When u is very small we again obtain 
the results of AE and T being small, this corresponds to a free particle; 
this is due to the fact that a free particle cannot absorb one photon 
and conserve momentum. 

There is also a very pronounced region of values U may have 
where we obtain quite large values for SE and [. This region corresponds 
to a loosely bound particle. This arises because the particle is not 
tightly bound and thus the electromagnetic field can exert more influence 


upon its energy. 


3. The Nuclear System. 


We take the normalized radial wave functions R, (r) and Re (4) 


to betl 
Rs)\F Gane. (8/3) 1/2 r exp(-ur?/2) 

(4-22) 
R(x) = a(v3 nyt exp Evr2/2) 


i.e. the 1S and 1P states for a harmonic oscillator potential well, 


eet 
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We notice that , andy have units of fm? if energy is in Mev and mass 
in Mev sec fmt, 
Substituting 4-22 into 4-8 we obtain 


orate 
SE rey age eat eka ec cae eo he ene 


The integral over r in 4-23 is given py*3 


* 33/2 (kr) 


Ee | exp{-(utv) 12/2} dr = (utv)7>/2 exp{-k?/2(ut+v). (4-24) 
0 (kr)3/2 


Substituting 4-24 into 4-23 we obtain for B(k) 


Peary 2 (y75)4/4 
m (uty) 


B(k) = - exp{-k2/2(utv)} . (4-25) 


Substituting 4-25 into 4-9 we obtain 


AE = - 16e* (v7: ee Pf k_exp{-k*/(utv)} dk | (4-26) 
am 0 (k+Ep- EQ) (utv) 2 


Substitution of 4-25 into 4-10 gives 


32e2 1/2 | exp{-(Ep-E,)*/(utv) } 


eS Fes (4+27) 
Cy ) ( R Pp hays 


Putting K and c into 4-26 and 4-27 where they should occur we obtain 


ie co a | lea (4-28) 
me? 0 [kt+(E,-Ep) /Uc] (uty) 
"a2 Hc (uty)? 
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then AE is given in units of Mev and Tf in units of sec’. 

We shall assume that ER = E, since the difference 
Ep = Ey ve. Since 4-28 and 4-29 are not symmetric we have two cases 
to consider. First, where y, the final energy, is fixed and y, the 


initial energy, varies; second where H is constant and v varies. 
de i fixed 


a. Weer G2) 


rT. axn(avyyon 


A v-1/yo/2 


2.  E£ixed 


a. Use 


r~ exp(-w fue! - 


AE % ie 


Nea pole 


Since both AE and f are finite everywhere we then obtain the following 


figures. 
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Figure 8. 


Figure 9. 
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Figure 10. 


Figure ll. 
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From figures 8-11 we notice that the dependence of AE and [I 
on the binding energy is of the same form as for the atomic case, section 
2. Thus, if the nuclear potential is very strong, i.e. large binding 
energies, AE and [ are quite small. This is due to the fact. that the 
particle is bound very tightly, hence the effect of the electromagnetic 
field, on the particle, is very slight. We notice that if the binding 
energy is very small, i.e. the particle is almost free, then AE and [ 
go to zero. This arises from the fact that a free particle cannot 
absorb one. photon and conserve momentum. 

There is also a range of binding potentials where we may have 
very large values of AE and f. This region corresponds to a loosely 
bound particle. The large values of AE and [ obtained in this region 
arises because the particle is influenced more by the electromagnetic 
field. 

While this analysis cannot discriminate between two similar 
potentials or give the exact shape of the potential it can be used to 


determine whether the potential is strong or weak. 
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CHAPTER V. SUMMARY AND CONCLUSIONS 


In this thesis we have considered the problem of a bound, 
charged particle interacting with an electromagnetic field. We have 
seen, 2-48 and 3-61, that the correction, to the frequency for the 
classical particle and the energy level for the quantum mechanical 
particle, are both of order e2. 

For the quantum mechanical model we obtained not only the 
level shift but also the line breadth, probability of transition, and 
the cross-section; these quantities all agreed, in the lowest order of 
e*, with previous well-known results of perturbation theory. This 
agreement with previous work indicates that the model used, 3-4, conforms 
somewhat to the actual physical world. 

In chapter IV we used the results obtained in chapter III to 
investigate the dependence of the line breadth and level shift on the 
binding potential. The results obtained from this investigation were 
in agreement with what one would expect, i.e. for a very tightly bound 
particle the line breadth and level shift were very small while for 
more loosely bound particles we saw that the line breadth and level 
shift became quite large. 

From the results of chapter IV we then can make some general 
statements about the nuclear potential. First, if the experimental 
values of line breadth and level shift of nuclear systems are small 
then the potential which binds the nucleons in the nucleus is very 
strong, somewhat like a very deep square well. And vice versa if the 


line breadth and level shift are large then the nuclear potential 
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would resemble a shallow square well. 


It is not clear, as yet, how much more information about the 


nuclear potential one can obtain by the type of analysis as done in 


chapter IV. However, since the atomic system is quite well understood 
this might serve as a guide for further and deeper studies into the 


connection between level shift, line breadth, and the binding potential. 
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